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Abstract 

It has been shown recently that a neutron placed in an external quasistatic electric field 
develops an induced electric dipole moment Pind due to quantum fluctuations in the QED 
vacuum. A feasible experiment which could detect such an effect is proposed and described 
here. It is shown that the peculiar angular dependence of pind on the orientation of the 
neutron spin leads to a characteristic asymmetry in polarized neutron scattering by heavy 
nuclei. This asymmetry can be of the order of 10^'^ for neutrons with epithermal energies. 
For thermalized neutrons from a hot moderator one still expects experimentally accessible 
values of the order of 10~^. The contribution of the induced effect to the neutron scattering 
length is expected to be only one order of magnitude smaller than that due to the neutron 
polarizability from its quark substructure. The experimental observation of this scattering 
asymmetry would be the first ever signal of nonlinearity in electrodynamics due to quantum 
fluctuations in the QED vacuum. 



PACS numbers: 12.20.Ds, ll.lO.-z, ll.lO.Lm, 29.27.Hj, 14.20.Dh 



1 Introduction 



Classical electrodynamics is well known to be a linear theory leading to the superposition princi- 
ple. At the quantum level the basic QED Lagrangian remains quadratic in the electromagnetic 
fields, so that the theory still appears to be linear. However, quantum fluctuations in the QED 
vacuum induce nonlinear effects that lead to a breakdown of the superposition principle [Ij. In 
particular, these QED fluctuations make the vacuum appear as an electrically and magnetically 
polarizable medium. The size of these corrections in nonlinear QED (NLQED) is very tiny, so 
that experiments with ultra-high intensity lasers have been proposed to search for these effects, 
e.g. e"'"e~ pair production from the vacuum [2]-l3], vacuum birefringence [l]-[5], light diffraction 
by a strong standing electromagnetic wave [i6] , and nonlinear Compton scattering [7] . A different 
proposal, involving quasistatic external electromagnetic flelds interacting with given electric or 
magnetic sources, has been made recently [8]-[9j. In |8] general expressions were obtained for 
the induced electric and magnetic flelds in such circumstances, and applied to the case of an 
electrically charged sphere in the presence of an external, quasistatic magnetic field. As a result 
of QED nonlinearity there appears an induced magnetic dipole moment, as well as corrections to 
the Coulomb field of the sphere. In spite of this being a dramatic effect, experimental detection 
appears very challenging. The complementary case of a purely magnetic dipole moment placed 
in an external, quasistatic electric field Eg was considered in [9]. The result is an induced electric 
dipole moment Pind, plus corrections to the magnetic field produced by the magnetic dipole. It 
was then suggested that the neutron could be used as a probe in the presence of large electric 
fields of order |Eo| — 10^'' V/m, such as present in certain crystals. A distinctive feature of this 
induced electric dipole moment, which should help in its detection, is its peculiar dependence on 
the angle between pind and Eq, or equivalently the angle between Pind and the neutron spin. 
In this paper we follow up on the experimental observability of such an induced electric dipole 
moment of the neutron. On the theoretical side we complete the analysis of [9| by computing 
the interaction Hamiltonian of the neutron immersed in a large external quasistatic electric field 
Eq, and an external, quasistatic, magnetic field Bq of ordinary strength. Given the nonlinearity 
of the problem one needs to check that (a) the magnetic interaction energy is of the usual form, 
(nonlinear magnetic corrections due to Bq are negligible), and (b) that the induced electric 
dipole does interact with the electric field Eq that generates it. The latter interaction energy is 
expected to have the standard functional form H^^t oc pind "Eq, albeit with an a-priori unknown 
coefficient which we determine. Next, we study the quantum behaviour of pind by means of the 
Heisenberg equation of motion. This is important for experiments based on potential changes in 
the Larmor frequency of the neutron spin around an external magnetic field due to the presence 
of Pind- We find no effect here, thus ruling out experiments of this type to detect an induced 
electric dipole moment of the neutron. Finally, we discuss in some detail a different approach 
based on neutron-nucleus scattering and conclude that this experiment offers an excellent op- 
portunity to observe such an effect. This is due to the peculiar angular dependence of pind- 
We find that for sufficiently large momentum transfers, a scattering asymmetry is induced with 
such particular characteristics that it would be easy to distinguish from other standard effects. 
The experimental discovery of such an asymmetry would be the first ever signal of a nonlinear 
effect in electrodynamics due to quantum fluctuations in the QED vacuum. 
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2 Induced electric dipole moment of the neutron 



An appropriate framework to discuss nonlinear effects induced by quantum fluctuations in the 
QED vacuum is that of the Euler-Heisenberg Lagrangian [T]. This is obtained from the weak 
field asymptotic expansion of the QED effective action at one loop order leading to 

4i = C(47-^ + 7^2) + ..., (1) 
where the omitted terms are of higher order in the expansion parameter In SI units 

C = — ~ 1.3 X 10-^2 (2) 

with OEM = e^/ (47reo/ic) the electromagnetic fine structure constant, nie and e are the mass and 
charge of the electron, respectively, and c the speed of light. The invariants and G are defined 
as 

j-=i(E2_c2B2)=-iF^,F'^^ (3) 

g = cB B = -\f^uF^\ (4) 

with F^y = d^Ay - dyA^ and F^"' = le^^P^Fp^. The so-cahed critical field E^, which plays the 
role of a reference field strength for the onset of nonlinearity, is given by 

he m 

This estimate is obtained by computing the electric field needed to produce an electron-positron 
pair in a spatial length of one Compton wavelength. For fields stronger than E^ the weak-field 
asymptotic expansion leading to Eq. (JT|) breaks down. In [8j general expressions were obtained 
for electric and magnetic fields induced by nonlinearity, to leading order in in the presence of 
external quasistatic weak fields (smaller than E^) and arbitrary sources. These induced fields 
are 

8 (x) = -^Vx / ^^Vy • f 47-mDm + -^mHm) , (6) 
B (x) = X / ^^Vy X ( - 4J-mHm + 7c^mDm ) , (7) 



27re§c2 7 |x-yl 

where Dm and Hm are the Maxwell (classical) fields produced by the arbitrary sources. Notice 
that these fields vanish as ^ — > (C — > 0). In the case of a current density uniformly distributed 
on the surface of a sphere of radius a, or equivalently, for a uniformly magnetized sphere of the 
same radius, the Maxwell, magnetic dipole type field, is given by 

^0 /3(m-er)er -m 2m \ 

Bd=4^| ^{^-o^ + ^Qia-r)"^, (8) 

where m is identified with the magnetic dipole moment of the source, and e.^ is a unit vector 
in the radial direction. Since the central expressions, Eqs.® and ([7]), were derived assuming 
E = cB < Ec the following constraint follows 



|m| 

he Ho 



3 < ^77^- (9) 
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For instance, if |m| = 0.96 x 10^^^ A m^, as for the neutron, then it follows that a > 10 fm. 
If this magnetic source is placed in an external, constant electric field Eq it has been shown [9] 
that there is an induced electric field of the dipole type 



f(x) 



1 p(V') 



IND 



X 



+ 



X 



(10) 



where "0 is the angle between the external electric field, lying in the x-z plane, and the magnetic 
dipole moment pointing along the z axis, i.e. Eg = |Eo| (sin^e^; + cos^l^ez), and m = |m|e2. 
The induced electric dipole moment p(?/')ind is given by 



p(V') 



IND 



CjUQ |m| |Eo 
lOvreoa^ 



36- 



Eq 

En 



49 



Eq 
En 



This induced electric field is of the electric dipole type in its radial l/|xp dependence, but it has 
a manifestly peculiar angular dependence. For instance, along the z axis, and unlike a standard 
electric dipole field, it has a non-zero component along bq that depends on the azimuthal angle 
(j). It also has a non-zero component along the direction of e^, as may be appreciated by writing 
the induced electric field in spherical coordinates (r, 9, 4>), i.e. 



5(x) 



C^o|m|2[E 



407r2e2a3 



— ^ I 2 [36 cos 9 cos ^ — 13 sin 9 cos (p sin ip] e^. 
x| 



+ [13 cos 9 cos (f) sin ^ + 36 sin 9 cos ^] eg — 13 sin (p sin ^/^ | 



:i2) 



In addition to the induced electric field Eq. (llOp . there is an induced magnetic field (a correction to 
the field produced by the magnetic dipole source), which can be derived from a vector potential, 
i.e. 13{x) = V X ^(x), where after a lengthy calculation one finds 



^(x) 



Cm 



^ I 4 |EoI^ (e^ X m) - 7 [m • Eo + 3 (Eq • e^) (m • e^)] (e^ x Eq) 
x| 



4vreo 

+7 (Eo • e^) (m x Eq) } 



(13) 



l + O 



fJ-o |m|- 
a^eo |Ec 



+ 



X 



Notice that while £ grows linearly with |Eo|, B depends quadratically on |Eo|. 
We proceed to discuss the interaction energy of the magnetic dipole source and its induced 
electric dipole with the external constant field Eq, and with an external uniform magnetic field 
Bq weak enough not to induce nonlinear effects, i.e. c|Bo| ^ Ec- Given the nonlinearity of the 
problem it is important to verify that the magnetic interaction Hamiltonian has the expected 
form — m • Bq, given the strength of Bq. In addition, the electric interaction energy of the 
induced electric dipole and the external field Eq is a-priori unknown. This need not be exactly 
of the form 

i^int = -^P-Eo, (14) 



as one would obtain for a linearly polarizable particle immersed in an external electric field, e.g. 
for a polarizable neutron on account of its quark substructure. In fact, the electric interaction 
Hamiltonian due to nonlinearity lacks the factor 1/2 as shown next. The canonical energy- 
momentum tensor is defined as 



dC 



tot 



Aot 5^., 



(15) 
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where the total Lagrangian density is Aot = — iii^^ ■, with L = e^F + C-^^i 
Eq.([T]). This equation can be rewritten as 



and /Igy given in 



(16) 



Since the last term on the right hand side above is the total divergence of an anti-symmetric 
tensor, employing the equations of motion 



one can define another energy-momentum tensor as 



(17) 



(18) 



Notice that this tensor is symmetric and gauge invariant except for the last two terms. The 
total energy density of the system is defined as the component 6^q, 



H 



tot 



dC 

— •E-/:-jA = DE-£-jA. 



(19) 



In general, for a given configuration of the fields the interaction Hamiltonian is defined as the 
difference of the total Hamiltonian with and without the sources. In a quantum theory it is 
defined as the difference of the total Hamiltonian evaluated at the fields in the interaction 
picture, with and without the external sources. Then, the interaction Hamiltonian i^mt, i.e. the 
volume integral of the interaction Hamiltonian density T-Lint is 



ifint = y"?^int d\ = - • Ad\ = - • (Ao + A) dV, 



(20) 



where A = Aq + A, with A given in Eq. (jl3p . and Aq is the vector potential associated with 
Bq, i.e. Bq = V X Aq, and Aq = ^ Bq x r. The current j corresponding to the magnetized 
sphere producing the field, Eq.([8]), is j = m x 5{r — a). In Eq. (pO]) the self energy of the 
magnetized sphere, independent of the external field, has been omitted. After performing the 
integration in Eq. (l20p one finds 



-m • Bo-p (V')iND -Eo, 



(21) 



which has the correct magnetic interaction term as in the linear theory. The electric interaction 
energy is of the expected form, but it involves a coefficient different from the case of linear QED 
as a result of nonlinearity. In the absence of the external magnetic field Bq, and using Eq. fllip . 
the interaction Hamiltonian becomes 



int 



C^o [mp [EqI^ 
lOvreoa-^ 



^36 - 49 sin^ ip) 



CfiQ |mp |Eo[^ 
lOTreoffl'^ 



'13 - 49cos2 V) 



(22) 



Notice the dependence of //int on a~^. It should be pointed out that in an experimental situation 
one would typically be interested in a point magnetic dipole. This source would produce very 
strong fields in its proximity so that the limit a — ?> would obviously not be allowed. Instead, 
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we assume that even in such a case the large distance solution for the fields is well described 
by the first order approximation to the effective Lagrangian C'^l in Eq.©. We also assume 
that this solution is robust against short distance modifications of the source as long as its 
symmetry is preserved. In this sense the parameter a is to be considered as a measure of our 
ignorance about the higher order corrections to this effective Lagrangian, something necessary 
when dealing with strong fields. The specific value of a will be discussed later in Section 4. The 
fact that H[^t depends on the orientation of m with respect to Eq through the angle ip can be 
used as a distinctive feature in the design of an experimental asymmetry as described below in 
Section 5. 

We consider next the quantum behaviour of pind using the Heisenberg equation of motion. To 
this end we consider a particle with magnetic dipole moment m related to the spin through the 
standard relation m = ghS, where g is the gyromagnetic ratio. Assuming that the dynamics 
of this particle is described by the Hamiltonian Eq. (j2T]) . and given that Pind oc |mp, iJjnt to 
first order in contains only quadratic terms in the spin, whose components are the dynamical 
variables of the problem. The effective Hamiltonian involving these dynamical variables must 
be symmetrized in order to ensure Hermiticity. Hence, the quadratic terms in the spin entering 
the Heisenberg equation of motion lead to the commutator 

[{Si, Sj} , Sk] = icjki {Si, Si} + ieiki {Sj,Si} . (23) 

For a spin 1/2 particle, such as the neutron, we have Si = ^0",, and {Si, Si} = ^Sn. In this case, 

1 I i 

[{Si, Sj} , Sk] = iejki-Sii + ieiki-5ji = - {ejki + (-ikj) = 0. (24) 

Therefore, dS/dt = so that if one is interested in the time evolution of a spin 1/2 particle, 
and Eq. (j2ip describes its effective Hamiltonian, we find no contribution from this leading order 
nonlinear correction. In other words, the precession of the spin is not affected. This is not the 
case, though, for spin-one particles. This unfortunate feature rules out experiments to detect the 
induced electric dipole moment of the neutron based on Larmor frequency changes. A different 
approach involving neutron scattering off nuclei is discussed next. 



3 Neutron-atom scattering amplitude and cross section 



Scattering of slow neutrons by a free atom can be described by a scattering amplitude in the 
Born-approximation, which in the center of mass system is given by (see e.g. [10| ) 



M f o 
/ (q, s) = / exp (iq • r) i^int (q, s) d r. (25) 

where M is the reduced mass 

M = (26) 
mn + niA 

with mn the neutron mass, and rriA the mass of the atom. The three-momentum transfer q is 
q = k — k', with k and k' the neutron wave vectors before and after scattering, respectively, and 
s is the neutron spin in units of h. The magnitude of q will be denoted as |q| = q in the sequel. 
The total Hamiltonian i^int involves all known interactions between the neutron and the atom, 
to which we add now the new interaction due to NLQED given in Eq. (j22p . Correspondingly, 
the total scattering amplitude can be written as 

/ (q, s) = /n (q, s) + /mag (q, s) + /c (q) + /pol (q) + /so (q, s) + /pv (q, s) + /ind (q, s) , (27) 



5 



where the various contributions are as follows. The term /n (q, s) is due to the hadronic in- 
teraction of the neutron with the nucleus, and is usually the dominant term. The amplitude 
/mag (Q) s) corresponds to the interaction of the neutron magnetic moment with the atomic mag- 
netic field (for atoms with unpaired electrons). This term can be of a similar size as /n (q, s). 
The next three terms arise from various electromagnetic interactions, i.e. /e (q) is due to scat- 
tering of the neutron charge radius by the electric charges in the atom, /poL (q) arises from the 
electric polarizability of the neutron due to its quark substructure, and /so (q, s) corresponds to 
the spin-orbit interaction of the neutron in the electric field of the nucleus. The term /pv (q, s) 
is a weak interaction, parity- violating amplitude which we list separately from /n (q, s) as it has 
a different dependence on neutron spin. Finally, /ind is the new component due to the induced 
electric dipole moment of the neutron, which we wish to isolate experimentally. 
The scattering amplitude, Eq. (j27p . enters the differential cross section for elastic neutron scat- 
tering by a single atom in the ground state, 

^(q,P) = (|/(q,s)|2), (28) 

which includes an ensemble average over nuclear and electronic spin degrees of freedom (if 
present), and the neutron spin. The incident neutrons are characterized by a polarization 
defined as P = 2 (s). In the absence of nuclear and electronic polarization of the atom, the case 
of interest here, the kinematic scattering variables are q and P. Experimentally, one determines 
neutron scattering cross sections using a sample containing a macroscopic number of atoms. 
Considering a single atomic species, the ensemble average in Eq.(|28p still has to account for the 
isotopic composition and the different states of total spin of a neutron scattering off a nucleus 
with non-zero spin. For slow neutrons with wavelengths much larger than the nuclear radius 
i?N the hadronic amplitude /n is practically independent of q. This is in the absence of nuclear 
resonances for thermal and epithermal neutrons, i.e. for the energy range of interest here. 
Scattering thus proceeds in an s-wave and is isotropic in the center of mass system. One defines 
a neutron scattering length operator as 

ON (s) = - lim /n (q, s) . (29) 

For a nucleus with spin hi one has 

(J + l)a++Ja, , 2(a+-a_) 
an (s) = + 21 + 1 ^ • I' (30) 

where a+ and a_ are the eigenvalues of ajsf (s) for the two states of total spin / it 1/2 (see e.g. 
|10]). For a sample with all nuclear species unpolarized, scattering by the ith isotope enters with 
statistical weight factors Wi+ = {li + 1) / (2/j + 1) and wi- = h/ (2/i + 1). The leading term in 
the cross section is then given by 

'^Ci\wi^\aij^\^ + Wi_\ai_\^ , (31) 



I |2 
ON 



where the bar indicates the averaging over isotopes and spin states, and Cj stands for the relative 
abundance of the ith isotope. In next-to-leading order the cross section contains interference 
terms between small amplitudes like /ind and a usually dominant coherent nuclear scattering 
length On, which for unpolarized nuclei is given by 



ON 



Ci [wi+ai+ + Wi-aiJ\ . (32) 
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Most scattering lengths on are found to be positive with typical values of a few fm. Neutron op- 
tical measurements determine a coherent bound scattering length b, related to the corresponding 
scattering length o of a free atom through 

\ + ^). (33) 

This relation includes the contributions — linig^^o /« due to all amplitudes (i = N, POL...) ap- 
pearing in Eq. (j27p . Lacking sufficiently accurate theoretical predictions for the nuclear part, 
however, one cannot even extract from b the sum of all non-hadronic components, which nor- 
mally contribute less than 1%. Instead, one needs to perform measurements for q ^ 0- Values 

— —2 

for b and the total neutron scattering cross section of an atom fixed in space, iTs,b = 47r6 , can be 
found e.g. in [11]. For later use we quote the values for lead with natural isotopic abundances 

b = (9.401 ± 0.002) fm, (34) 
(7s,b = (11.187 ± 0.007) X lO^^^'cm^. 

For low neutron energies one also has to take into account interference effects of the neutron 
waves scattered from different atoms, as e.g. in condensed-matter studies. Classical examples are 
Bragg scattering by single crystals and measurements of phonon dispersion relations. However, 
for sufficiently large momentum transfer q as considered here, interatomic interferences and 
eventual nuclear spin correlations between different atoms can be neglected. We thus consider 
the cross section in the center of mass system as given by 

^ ^ ^ + (l/so (q,s)p) + ... - ^ (2Re [Wfj (q,s)]) . (35) 

i 

where the dots stand for the remaining contributions of squared amplitudes from Eq. ()27p . and 
the sum is over j = e, MAG, POL, SO, PV and IND. The formulas to transform this expression 
to the laboratory reference frame can be found in Ref. |10j . As long as the nucleus is free 
to recoil, the total cross section does not change and changes in the angular distribution of 
the scattered neutrons appear only at order m^/rriA- Since we are not interested in angular 
distributions and will only consider atoms much heavier than the neutron, we may use the 
scattering cross section as given above in Eq. (l35p . 



4 Scattering amplitude due to the NLQED induced electric 
dipole moment 

We now turn to the calculation of the new amplitude /ind due to the NLQED induced electric 
dipole moment Pind given in Eq. (|lip . The magnetic moment of the neutron, m, can be written 
as 

m = /XnfT, (36) 

where 

/Xn = -9.662 X 10^27 A m^ , (37) 
and o" = 2s are the Pauli matrices. From Eq. Q there follows the lower bound 

a>7.6fm. (38) 
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According to Eq. (|2ip . pind interacts with the atomic electrostatic field, which for simplicity we 
consider as given by a pointlike nucleus with electric charge Ze, 



En 



1 Ze 
47reo 



(39) 



As discussed later, one can neglect electric field shielding due to the atomic electrons. Using 
Eqs.dSSD, (122]) and (gSD one obtains 



(40) 



/iND (9, a, /?) = g^o^^^^y [-13 h (g, R) + 49 h (q, R, /3)] 



where /3 is the angle between s and q. The two integrals Ii {q, R) and I2 {q, R, (3) can be easily 
calculated analytically in polar coordinates with q along the polar axis. They are given by 



and 



h{q.R,P) 



h{q,R) = 
exp (iq • r) 



exp (iq • r) 3 
d r. 



(cos /3 cos 9 + sin /3 sin 6 cos ip)'^ d^r. 



(41) 



(42) 



The radial integration extends over all space, excluding a sphere of radius R around the nucleus. 
For a heavy nucleus like lead, electric fields as strong as 10^^ Vm~^ exist close to the nuclear 
surface. This exceeds by far the critical field, Eq.([S]), beyond which higher order terms in the 
one-loop effective Lagrangian in Eq.([T|) become important [12J and thus cannot be neglected. 
Therefore, R has to be much larger than the nuclear radius Rt<{, and we choose it here as the 
distance from the nucleus to where the critical field is reached, i.e. 



Er 



1 Ze 
4^^ 



(43) 



For lead isotopes with Z = 82 one has R ~ 300 fm. The integrals Eqs.(41) and (42) can be 
solved analytically with the result 



2lT 



h{q,R) = — {cos {qR) + 



sin (qR) 
qR 



+ 



vr 



Si (qR) -- 



qR 



(44) 



An 

'r 



where Si (x) is the Sine integral, and 



l2iq,R,P) 



vr r 1 

4^? l(^ 



-qR+-iqRf 



2 + 3 {qRf + ( 6 + {qRf ) cos (2/3) cos (qR) 



2 - 3 {qRf + (6 - {qRfj cos (2/3)J sin (qR) 
-(7i?(3 + cos(2/3))(|-Si(gi?))} 



(45) 



47r 

m 



— (3 + cos (2/3)) qR+ — {2 + cos (2/3)) {qRf 
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Using Eqs. (|44p and ()45p in Eq. (j40p one finds the final expression for the scattering amphtude 
due to the NLQED induced electric dipole moment 



/iND (q, R, a, 13) 



CM/ioK^^e^ r 49 [1 + 3cos (2/3)] [qRcos {qR) - sin (qR)] 



320TT^h^ela^R 
+^[43 + 49 cos (2/3)] 



cos (qR) + 



2{qR)' 
sin (qR) 
qR 



(46) 



+ (Si (qR) 



qR 



2ATT'^h'^ela?R 



— [43 + 49 COS (2/3)] qR 



+100 



49 COS (2/3)] (gi?)^ 



The scattering amplitude /ind exhibits a welcome peculiar dependence on the angle /3 between 
the neutron spin s and the three-momentum transfer q. This dependence introduces an asymme- 
try which for suitable experimental conditions is essentially free from background contributions 
due to other well known effects. The largest effect is obtained by evaluating /ind at /3 = 
and at /3 = 7r/2. This feature will then play a crucial role in the experimental detection of the 
NLQED induced electric dipole moment of the neutron, as discussed in the following section. 



5 Scattering asymmetry due to the NLQED induced electric 
dipole moment 

We define the scattering asymmetry as 

(da/dO)|| - (da/dl))^ 

where the differential cross section, Eq. (j35p . is evaluated for two neutron polarization states Py 
and Pj_, parallel and perpendicular to the scattering vector q, respectively. The interference 
term between the coherent nuclear amplitude and the amplitude of interest leads to 

47r 

^IND (g, -R, a) = W (/lND± - /iNDIl) P-, (48) 

where /ind|| = /ind (g, ^, a, 0) and /ind± = /ind (g, ^, a, vr/2), P = \P\\\ = |P_l|, and = 
<^s,hM'^ /m'^ is the total scattering cross section of the free atom. We argue below that in a well- 
designed experiment possible influences of interference terms other than between on and /ind 
are negligible. Hence, A{q,R,a) ~ Ansfj:, {q, R,a), so that the asymmetry defined in Eq . ()^7|) 
should allow for a detection and determination of the new amplitude. Using the values for 
natural lead from Eq. (l34p in Eq. ()48p one obtains 



4 / n N /lND± — /iNDII „ . .„N 

Aiq,R,a)^^^^^R (49) 

Prom Eqs.dMI) and (jH]) it follows that A{q,R,a) oc xil^) /R-i where X is a function of the 
dimensionless parameter qR. The maximum of A (g, R, a) occurs for (see Fig. 1) 

qR = 1.68. (50) 
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Figure 1: The scattering asymmetry A{q, R = 300 fm, a = 7.6 fm), Eq. (j49p . normalized to the 
neutron polarization P as a function of the dimensionless variable qR. 

Using the value of a given in Eg . (1381) the asymmetry becomes 

A {5.6 X 10^^m~\ 300 fm, 7.6 fm) = 1.4 x lO^^P, (51) 

a result which appears experimentally accessible. Neglecting nuclear recoil, valid in good ap- 
proximation for neutron scattering off a heavy target, one may use the relation 

q = 2k sin (52) 

where is the angle between k and k', and /c = |k| = 2.197 x 10~^ fm~^y^ E{eV) , with E 
the neutron kinetic energy in eV. In a backscattering geometry, i.e. for ~ tt, the maximum 
asymmetry, Eq. ()5ip . is obtained with epithermal neutrons of energy E ~ 165 eV. The result for 
the asymmetry depends explicitly on the radii R and a of spheres centered on the nucleus and 
on the neutron, respectively. These parameters play the role of separating the short-distance 
physics close to the electromagnetic sources from the long-distance effects involving fields weak 
enough for the Euler-Heisenberg approximation to be valid. In this regard it is important 
to point out that electric and magnetic fields induced by quantum fluctuations in the QED 
vacuum involve various multipolarities [8]. For instance, in the case of a neutron in an external, 
quasistatic electric field |Eo| < Ec, the induced electric field Eq. p^ has a dipole type term 
of order 0(|x|~^), as well as a higher multipole of order Odxp^), while the induced magnetic 
field involves terms of order C(|x|~^), 0(|xp"'), and 0(|x|~^), with x = r — rn- These higher 
order multipoles can be safely neglected in the interaction energy as long as the weak-field 
approximation remains valid. However, it is not clear what happens at distances closer to the 
nucleus or to the neutron. This problem is similar to that of the separation into far and near field 
regions around a localized charge/current distribution in classical electrodynamics, where the 
lowest-order multipole provides the long-distance solution. Although one cannot compute the 
interaction energy due to QED vacuum effects stemming from the regions r < R and |r — rn| < a, 
their contribution to the asymmetry would presumably have a different dependence on s and q. 
For large q it might smear out the oscillations appearing in Fig. [H For small q , corresponding 
to small spatial resolution in probing the QED vacuum, the asymmetry should not be affected 
much by the short-distance physics. This statement is underlined by the fact that at leading 
order in qR the /3-dependent term in /ind, Eq. ([l6]) . does not depend on R. For qR <^ 1, one 
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thus obtain a prediction for the asymmetry which should be robust against variations in the 
choice of R, i.e. 

A{q«R ,R,a)^Aiq,a) = —- ^^,^3^3 qP . (53) 

For a neutron energy of 1 eV in a backscattering geometry one expects 

A (4.4 X 10^^m"\ 7.6 fm) = 2.2 x 10~^P . (54) 

Prom a practical point of view the polarization of epithermal neutrons with more than 100 eV 
requires a spin filter of polarized protons. This is technically demanding if one wishes to polarize 
a beam with a diameter of several cm. The measured energy-dependent neutron polarization 
cross section of a polarized-proton spin filter for the energy range of interest may be found in Ref. 
|13] . On the other hand, neutrons with energies up to ~ 1 eV are available from a hot moderator 
in a reactor neutron source with much higher intensity than epithermal neutrons. They can be 
polarized using magnetic monochromator crystals, or by a spin filter of polarized ^He gas jl4j 
which has a polarization cross section proportional to For the fluxes available at the Institut 
Laue-Langevin (ILL) in Grenoble an asymmetry as in Eq.(|54p appears experimentally accessible 
within a few days of beam time. 

To conclude this section we stress that the neutron scattering asymmetry due to nonlinear QED 
has two characteristic properties which should make it rather easy to detect and distinguish from 
other effects. First, the asymmetry attains its maximum value for perpendicular orientations 
of the neutron polarization and vanishes for opposite orientations. This is in contrast to most 
ordinary asymmetries which become maximal for opposite orientations. Second, Aj^d (q, R, a) 
exhibits a characteristic q dependence with a broad maximum around the value of q given in 
Eq. ([50ll . These features are discussed in more detail in the sequel. 



6 Analysis of background asymmetries 

In this section we study the contributions to the asymmetry A (g, R, a), Eq. (j47l) . from the various 
ordinary scattering amplitudes defined in Eq. ()27p . The neutron spin-dependent amplitudes can 
be written as 

/(q,s) = /o(q) + /i(q)[s-w(q)] (55) 

where /o (q) is spin independent, and w is a vector not correlated with the neutron spin. In the 
case of the weak amplitude fpy the vector q must be replaced by k. For instance, for the nuclear 
amplitude in Eq. (j30p . w is independent of q and given by the nuclear spin I. It can be shown 
in general that the terms proportional to ((s • w)^) in the differential cross section, Eq. (l35p . are 
all independent of the neutron polarization and therefore cannot generate an asymmetry. In 
principle these terms influence the size of A {q, R, a) through the total scattering cross section 
(Ts, Eq. (jl8|) . For scattering angles G ^ the pure spin-orbit cross section, quadratic in the 
amplitude /sO) might become large enough to have an impact on as- However, for sufficiently 
large Q, and in the absence of nuclear and electronic polarization, corrections to Us due to 
squared-amplitude terms can be safely neglected. The interference terms between the nuclear 
and the other scattering amplitudes in Eq. ()35p may however affect the asymmetry A{q,R,a) 
through their dependence on the neutron polarization. This requires careful consideration, and 
we start with the amplitude /so for spin-orbit scattering. It originates in the interaction of the 
neutron magnetic moment with the magnetic field present in the neutron rest frame due to its 
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motion through the atomic electric fields. Its expression is (see e.g. [15j ) 



/so (q, s) = i— cot (G/2) ^eZ [I - F {q)] (s ■ n) , (56) 

where eZ [1 — F [q)] is the Fourier transform of the electric charge density of the atom. This 
term involves the nuclear charge Z and the atomic form factor F [q) normalized to F (0) = 1. 
This form factor is measured e.g. in X-ray scattering off atoms, and is a real function of the 
momentum. The unit vector n points along k x k', so that the amplitude can contribute only 
if the neutron polarization has a component out of the scattering plane. The asymmetry due to 
the spin-orbit interaction is given by 

Aso = cot (e/2) ^eZ [l-F (,)] (P,, • n - • rr) . (57) 

The imaginary part of on above is due to nuclear absorption and can be calculated using the 
optical theorem. Lead nuclei absorb neutrons only weakly so that Im on — (47r)~^ ka^ for 
neutron energies of interest here. For a scattering angle = 7r/2 a neutron kinetic energy of 
E ~ 330 eV would be required to observe the maximum asymmetry according to Eq. (l51jl . In 
this case the atomic form factor F {q) ~ 0, and a maximum asymmetry ^so should be observed 
for P_L perpendicular to the scattering plane, i.e. Pj_ • n = P (while Py • n = by definition). 
As a result, for E = 330 eV one has Aq,o ~ 4.8 x 10~^P, while for E = leV together with the 
conservative value F (q) = 0, the asymmetry becomes ^so — 2.6 x 10~^P. These values are 
already smaller than the corresponding values of Aind, but with suitable experimental settings 
they can be reduced even further. Notice that the condition P_l • q = required for /ind± can 
be realized for different orientations of P^, while also fulfilling Pj^ • n = 0. Hence, choosing a 
backscattering geometry, for which cot (0/2) <^ 1 together with P^ • n ~ 0, and using realistic 
assumptions about the experimental definition of the directions of q and P, one can easily sup- 
press Aso by a factor 50. Hence, the impact of Aq,o can be kept well under the 1% level. 
Next, the amplitude on of the neutron-nuclear interaction in Eq. (l30p . when squared, gives rise to 
interference between its spin-dependent and spin-independent parts. After ensemble averaging 
this becomes proportional to P and to the nuclear polarization P^. In thermal equilibrium, 
Pn = tanh(^Ni?/ (kBT)) for nuclei with magnetic moment /^n in a magnetic field B (/cb is the 
Boltzmann constant). If the target is at room temperature, and given that no magnetic field is 
needed at the position of the sample, the P-dependent cross section is orders of magnitude too 
small to have an impact on the asymmetry. 

We consider next the amplitude /maGj which is due to the interaction of the neutron mag- 
netic moment with the magnetic field produced by unpaired atomic electrons of paramagnetic 
contaminants in the sample. The operator structure of this amplitude is given by 

/mag oc s • [eg X M (q) X e,] , (58) 

where M (q) is the Fourier transform of the total (spin and orbital) magnetization of the atom, 
and Bq = q/q. The interference term with cin thus involves the neutron polarization and the 
sample-averaged magnetization. It would only influence the asymmetry if (1) the ratio B/T is 
sufficiently high to result in a sizable magnetization, (2) paramagnetic centers are sufficiently 
abundant, (3) the two neutron polarization states in the asymmetry have different projections 
perpendicular to the scattering plane (which is a consequence of the term in brackets in Eq. (j58p ). 
and (4) measurements are performed for sufficiently small q, where the magnetic form factor still 
has a sizable value. Regarding the latter, even for the smaller value of q envisaged in Eq. (j54p . 
the magnetic form factor leads to a strong suppression of /mag- Hence, with the conditions 
1, 2 and 3 under experimental control one can safely disregard magnetism as a source of an 
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asymmetry. 

Finally, the parity-violating amplitude /pv due to the hadronic weak interaction may lead to 
a different type of asymmetry which has indeed been observed in neutron transmission experi- 
ments. Effects depend on the neutron helicity, hence w = k in Eq. (l55p . with a complex coefficient 
to describe both parity violating spin rotation and transmission asymmetry. The amplitude is 
normally so small that it requires special efforts to detect it. For thermal neutrons, transmission 
asymmetries for longitudinally polarized neutrons [16] have typical sizes of a few times 10~^. 
However, for neutron energies in the vicinity of p-wave resonances of complex nuclei a strong 
enhancement due to the weak nuclear interactions may appear. A prominent example is the 
transmission asymmetry of 7% found at the p-wave resonance of 0.76 eV in ^'^^La [17J. However, 
no effect sufficiently strong to affect the asymmetry ^ind is known for lead in the relevant energy 
range. In addition, an experimental test can easily be performed. In fact, taking the neutron 
polarization P parallel and anti-parallel to q, i.e. f3 = and /? = vr, it follows from Eqs. (l46p 
and (|48p that ^ind = for these two polarization orientations. In contrast, for fpy one has 
^PV oc sin (0/2), which could be measured separately and corrected for if the need arises. 



7 Comparison of the NLQED amplitude with ordinary electric 
amplitudes 

The electric amplitudes /poL and /e do not generate any known scattering asymmetry. However, 
owing to their characteristic g-dependences a comparison with the NLQED amplitude /ind is 
needed. Like /inDj the amplitude /poL due to the electric polarizability of the neutron, «„, is 
induced by the nuclear electric field, Eq. ([39]) . In SI units is defined by p = 47reoanEo, so that 
its dimension is [an] = m^. The calculation of /poL follows from Eq. (j25l) with the interaction 
energy given by Eq. (jl4p . This leads to the integral Eq. (|4ip with the result given in Eq. (j44p . 
However, the lower limit of the radial integration is now different from that for /ind- In fact, 
this lower limit can now be extended down to the nuclear radius R^, since for r > the 
neutron probes only the long-range electric forces. For r < R^ the electric interaction is small 
in comparison with the nuclear force, so that in early calculations [20j-[21j it has simply been 
included in the nuclear amplitude. In SI units and for the electric field given in Eq. ()39p . the 
dependence of /poL on q is given by 

The term linear in q is characteristic of the dependence of the Hamiltonian, and it also 
enters the interference term in the cross section. This feature has been exploited in the past to 
measure «„• Conflicting results from experiments performed during more than three decades 
show that a proper assessment of all systematic errors has been difficult (see e.g. the table in 
Ref. [22]). The most recent result [23], derived from energy-dependent neutron transmission 
through a ^osp]^ target, and reporting the smallest uncertainty is 

On, exp = (1.20 ± 0.15 ± 0.20) X IQ-^ fm^ (60) 

Calculations using quark bag models [24J agree with this result. An early estimate of Breit 
and Rustgi [25\ using data on pion photoproduction already indicated that < 2 x 10~^fm^. 
These authors also analyzed other effects which might mimic a signal from the neutron electric 
polarizability. From an estimate of vacuum polarization effects close to the nucleus, and using 
the Uehling potential [26], they concluded that this contribution to neutron scattering can be 
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safely neglected. Turning now to the question of shielding of the nuclear charge by the atomic 
electrons, one notices that this might quench the amplitude /ind- Since /poL and /ind both 
depend quadratically on the electric field Eq, one can draw parallels with the analysis of /poL- 
We recall that in our calculation of /ind one needs to exclude a spherical region of radius R 
around the nucleus inside which the weak-field expansion of the Euler-Heisenberg Lagrangian 
breaks down. This procedure was followed (for different reasons) in the early calculations of /poL 
|20j-[21]. |25j where the nuclear region with radius R-^ was excluded. A more recent analysis 
which does not rely on a simple model for the nuclear charge distribution gives [2^ 



1 [SMZ^e^ , , 

/pOL(g^0) = - a/-T2 «n, (61) 

47reo \ IT n'^ 

where is the root mean square charge radius of the nucleus. This shows that slow neutrons 
are indeed insensitive to details at this length scale. Notice that this result is nearly identical to 
the leading term in Eq. (|59p after replacing R^ by r^. In the derivation of Eq. (j6ip the following 
intermediate result was obtained in [27] 

/pOL(g^0)oc/ \F^{K)-F{K)\^dK= |FN(K)|2d/t[l-0(i?N/i?A)], (62) 

JQ Jo 

where -Fn (f^) and F (k) are the charge form factors of the nucleus and of the electron distribution 
in the atom, respectively. With R^/Ra — 10~^, shielding of the nuclear charge can be neglected 
in /poL- Even in the limit g — )• the neutron feels the full unscreened nuclear charge as far 
as the electric polarizability is concerned. For the NLQED induced electric dipole moment it 
follows that, with R ~ 300 fm for lead, the corresponding correction term of order 0{R/R\) 
is much larger than the one of order O (Ri^/Ra). However, R/Ra ^ 10~^ is still small enough 
so that one can neglect shielding of the electric field in the region around the nucleus. We now 
compare the two amplitudes /poL and /ind in the limit g — 0, i.e. their respective contributions 
to the neutron scattering length. Setting i^N = 1.2fmA^^^ and using Eq. (160p one can estimate 
the leading order term in Eq. (j59p as 

/poL(g^O) ~0.04fm. (63) 

Prom Eq. (j46p the corresponding leading term in /ind is 

1 MZ^e^ ^ ^ 

/iND(g^0) = ^^^aiND, (64) 



where 

C/Uo/in 3.3 fm^ 
QTT^ela^ ^ (a [fm]) 

With the value of a from Eq.Q one obtains 



ajND = 7.5 X 10"^ fm, (66) 

which is larger than a^, Eq. ()60p . However, the impact of this polarizability on the amplitude 
/ind for q — > is suppressed with respect to /poL due to R^ i^N- In fact, for lead one obtains 

/ind (g ^ 0) - 0.006 fm. (67) 

This result, though, is not even an order of magnitude smaller than /poL (<? — ^ 0), Eq. ([63]) . It 
would thus contribute about 5 x 10~^^ eV to the neutron optical potential of solid lead, which 
is quite substantial given the high precision of some neutron optical methods. The contribution 
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Figure 2: Amplitude (/ind {q,R = 300 fm, a = 7.6 fm)) for unpolarized neutrons, Eq. ([69|) . as a 
function of the dimensionless variable qR. 



of the NLQED electric dipole moment to the total cross section is (Tind — — Svra/iND (<? — ^ 0), 
where on ^ a in Eq. ([35|) has been used. Numerically this becomes 

o-iND (g ^ 0) = -0.014 X 10~2^cm2. (68) 

We discuss next the amplitude /e which describes the interaction of the electric charges of the 
atom with the internal charge distribution of the neutron as characterized by its mean squared 
charge radius. The amplitude for a bound nucleus is given by /e = —beZ [1 — F (q)], with the 
atomic form factor F (q) given in Eqs. (j56p and (j62p . and the neutron-electron scattering length is 
be — —1.35 X 10"'^ fm as determined from measurements of the total cross sections of lead and bis- 
muth at different neutron energies |18j-[19j. For lead and sufficiently large q, fe = 0.11 fm, which 
leads to a contribution to the total scattering cross section ^27j Ue — —STrabeZ ~ 0.25 x 10~^^ cm^. 
It is also interesting to notice that since /e (^ — > 0) ^0 this amplitude does not contribute to 
the neutron scattering length. The atomic form factor changes significantly at small q where 
interference of neutron waves from different atoms cannot be neglected. Hence, the macroscopic 
state of the sample enters crucially in the analysis of scattering data. In contrast, in the case of 
/ind where larger effects show up at much higher values of q, interatomic interferences do not 
play any significant role. 

To conclude this section we discuss the contribution of /ind as a potential background in mea- 
surements of the amplitudes /poL and /e performed with unpolarized neutrons. For the polar- 
ization averaged /ind one obtains 



1 r 

(/ind {q, R,a)) = - /ind {q, R, a, P) sin (3 6.(3 
^ Jo 

(MfiofilZ^e^ ( sm{qR) 



To analyze the impact on /c one may approximate relativistic Hartree-Fock results for F (q) by 
the simple function 1 +3{q/qo) [27\. With F (qo) = 1/2 the momentum go provides a 

scale at which significant changes in /e take place. For lead, qo = 8.3 x 10^*^ m^^, and 

/e('Zo)-/c(0) = 50xlO-3fm. (70) 
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In contrast, the change of (/ind) is much smaher due to its milder g-dependence and its smaller 
magnitude (see Fig. [2]), 

(/ind (go)) - (/ind (0)) = -0.11 X 10-3 fm. (71) 

Since the precision of the best measurements of be is at the level of a few percent, potential 
background due to (/ind) is negligible. A similar argument leads to the same conclusion for the 
determination of from /poL- 



8 Conclusions 



Many, if not most proposals to detect nonlinear effects due to quantum fluctuations in the QED 
vacuum rely on experiments involving lasers of ultra-high intensities [2J-[7J. These intensities, 
though, are at least two orders of magnitude below current values. An alternative approach has 
been discussed in this paper, based on the theoretical prediction of an induced electric dipole 
moment of the neutron, PinDi in an external quasistatic electric field [9]. The peculiar features 
of this dipole moment, particularly its dependence on the angle between Pind and the neu- 
tron spin, suggests the definition of an asymmetry which could be detected in the scattering of 
polarized neutrons from heavy nuclei. We have introduced this asymmetry and discussed all 
possible sources of background asymmetries. We have also compared the new NLQED ampli- 
tude with ordinary electric scattering amplitudes, particularly the one due to the polarization of 
the neutron in an electric field due to its quark substructure. The conclusion from this detailed 
analysis is that the asymmetry due to NLQED should be observable using epithermal neutrons, 
and even using thermalized neutrons from a hot moderator. This would be the first ever ex- 
perimental confirmation of nonlinearity in electrodynamics due to QED vacuum fluctuations. 
The numerical predictions for the asymmetry made in this paper were calculated using definite 
values for the parameters R and o. These were derived from the condition that the electric and 
magnetic fields should be below their critical values, beyond which the weak-field expansion of 
the effective Lagrangian breaks down. While the value of the asymmetry A for small q does not 
depend on R, it does depend on a as seen from Eq. (j22p Hence, the numerical results given here 
should be correct up to a numerical factor of order one. 
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